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Multiloop calculations in HQET 

A.G. Grozin 3, 

a Budker Institute of Nuclear Physics, Novosibirsk 

Recently, algorithms for calculation of 3-loop propagator diagrams in HQET and on-shell QCD with a heavy 
quark have been constructed and implemented. These algorithms (based on integration by parts recurrence 
relations) reduce an arbitrary diagram to a combination of a finite number of basis integrals. Here I discuss 
various ways to calculate non-trivial bases integrals, either exactly or as expansions in e. Some integrals of these 
two classes are related to each other by inversion, which provides a useful cross-check. 


I presented a review talk about multiloop cal¬ 
culations in HQET at this conference in Pisa 
in 1995 0]. Methods of calculation of two- 
loop propagator diagrams in HQET [p| and on- 
shell massive QCD , based on integration by 
parts ||, were discussed there. Recently, three- 
loop HQET |J and on-shell 0 algorithms have 
been constructed. Here I discuss this substantial 
progress. 


1. Three-loop massless diagrams 


First, I briefly remind you the classic method 
of calculation of 3-loop massless propagator dia¬ 
grams. There are 3 generic topologies of such di¬ 
agrams. They can be reduces, using integration 
by parts, to 6 basis integrals ||. This algorithm 
is implemented in the package Mincer Q (first 
written in SCHOONSCHIP || and later rewrit¬ 
ten in FORM [|I(| ), and in the package Sheer [ [U] 
written in REDUCE (T2)JT^. 

Four basis integrals are trivial. One is a two- 
loop diagram with a non-integer power of the 
middle line. It can be found as a particular case 
of a more general expression |lj,|l^| for the the 
two-loop diagram with three non-integer powers 
via a hypergeometric 3 F 2 function of the unit 
argument, with indices tending to integers at 
e —> 0. There is a rather straightforward algo¬ 
rithm for expanding such functions in £, with 
coefficients expressed via multiple ((-values. I 
have implemented it in REDUCE in the sum¬ 
mer of 2000, some results produced by this pro¬ 


gram are published in [[h| . It is clearly presented 
as Algorithm A in jlTj]; this paper also contains 
other, more complicated, algorithms. The algo¬ 
rithms of 0 are implemented in the C-)—b li¬ 
brary nestedsums [Q based on the computer- 
algebra library GiNaC [|l9|. This implementa¬ 
tion is very convenient; unfortunately, it requires 
one to install an outdated version of GiNaC. The 
Algorithm A seems to be also implemented in 
FORM pdf , but I could not understand how to 
use it. Using my REDUCE procedure or nest¬ 
edsums |l8| , one can quickly find as many terms 
of expansion of this basis integral in e as needed, 
in terms of multiple ((-values. They can be ex¬ 
pressed, up to weight 9, via a minimum set of 
independent ((-values, using the results of p2| , p0[ . 
The two-loop diagram with a non-integer power 
of the middle line can also be expressed ]2l[ via an 
3 F 2 function of the argument —1. Expanding this 
expression in £ (say, using nestedsums 0), we 
encounter more general Euler-Zagier sums, which 
were also considered in 221. Reducing them to 
the minimal basis, we obtain, of course, the same 
£-expansion of our basis integral. 

Using this expansion and integration-by-parts 
relations, it is easy to recover the well-known re¬ 
sult for the 3-loop ladder diagram, which is finite 
£ = 0: 20((5) + O(s). The last and most difficult 
basis diagram is non-planar. It is also finite at 
£ = 0. Using gluing of its external vertices ||, 
one can easily understand that it has the same 
value 20((5) at £ = 0. There is no easy way to 
find further terms of its £- expansion. 
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Figure 1. Inversion relations 


2. Three-loop HQET diagrams 

There are 10 generic topologies of 3-loop 
HQET propagator diagrams. They can be 
reduced, using integration by parts, to 8 ba¬ 
sis integrals Q. This algorithm is imple¬ 
mented in the RE DUCE package Grinder p|], 
available at http://www-ttp.physik.uni- 


karlsruhe.de/Progdata/ttp00/ttp00-01/. Five 


basis integrals are trivial. Two can be expressed 
via 3 F 2 hypergeometric functions of the unit 
argument [|23|,p| . Their expansions in e can be 
obtained in the same way as in the massless case, 
the results are presented in |h]]. The last and 
most difficult basis integral was found in [|o) up 
to the finite term in e, using direct integration 
in the coordinate space. More terms of its e- 


the FORM package SHELL2 Q. Two basis in¬ 
tegrals are trivial, and the third one is expressed 
via two 3 F 2 hypergeometric functions of the unit 
argument. However, some of their indices tend 
to half-integers at s —► 0 , and the algorithm of 
expansion in e discussed in Sect. |l] is not appli¬ 
cable. This approach was used for QCD/HQET 
matching of heavy-light quark currents |25| and 
chromomagnetic interaction |p(i| . 

The case when there is another non-zero mass 
was systematically studied in j2?j]. There are 
4 basis integrals, 2 of them trivial, and 2 are 
expressed via 3 F 2 hypergeometric functions of 
the mass ratio squared. Finite parts at e —> 0 
are expressed via dilogarithms. More terms of 
expansions of the general results p7| in e were re- 


expansion were recently obtained in ||24|| using 

is available at http://wwwthep.physik.uni- 

inversion, as explained in the next Section. 

mainz.de/Publications/progdata/mzth9838/ 


3. Three-loop on-shell diagrams 

Calculations of on-shell diagrams with massive 
quarks in QCD are necessary for obtaining coeffi¬ 
cients in the HQET Lagrangian and 1/m HQET 
expansions of QCD operators by matching. There 
are 2 generic topologies of 2 -loop on-shell propa¬ 
gator diagrams with a single non-zero mass. They 
can be reduced, using integration by parts, to 3 
basis integrals. This algorithm is implemented 
in the REDUCE package RECURSOR @ and 


Mm.red. 

There are 11 generic topologies of 3-loop on- 
shell propagator diagrams with a single non-zero 
mass (10 of them are the same as in HQET, and 
one involves a heavy-quark loop). They can be 
reduced, using integration by parts, to 18 basis 
integrals ]?j]. This algorithm is implemented as 
the FORM package SHELL3 (?[]. The basis inte¬ 
grals are mostly known from QED [p9| . 

Some on-shell diagrams are related to HQET 
ones by inversion of Euclidean integration mo¬ 
menta. One- and two-loop relations were pre- 
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sented in Three-loop relations are shown in 
Fig. U The second of them was used in Q to 
relate the the convergent ladder HQET diagram 
at e = 0 to the known on-shell ladder diagram. 
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